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Abstract 

The Georgia Tech/Goddard Global Ozone Chemistry Aerosol Radiation and Transport 
(GOCART) model is used to simulate the atmospheric sulfur cycle. The model uses the 
assimilated meteorological data from the Goddard Earth Observing System Data Assimilation 
System (GEOS DAS). Global sulfur budgets from a 6-year simulation for SO 2 , sulfate, 
dimethylsulfide (DMS), and methanesulfonic acid (MSA) are presented in this paper. In a 
normal year without major volcanic perturbations, about 20% of the sulfate precursor emission 
is from natural sources (biogenic and volcanic) and 80% is anthropogenic: the same sources 
contribute 33% and 67% respectively to the total sulfate burden. A sulfate production efficiency 
of 0.41-0.42 is estimated in the model, an efficiency which is defined as a ratio of the amount of 
sulfate produced to the total amount of S0 2 emitted and produced in the atmosphere. This 
value indicates that less than half of the S0 2 entering the atmosphere contributes to the sulfate 
production, the rest being removed by dry and wet depositions. In a simulation for 1990, we 
estimate a total sulfate production of 39 Tg S yr -1 , with 36% and 64% respectively from in-air 
and in-cloud oxidation of S0 2 . We also demonstrate that major volcanic eruptions, such as the 
Mt. Pinatubo eruption in 1991, can significantly change the sulfate formation pathways, 
distributions, abundance, and lifetime. Comparison with other models shows that the 
parameterizations for wet removal or wet production of sulfate are the most critical factors in 
determining the burdens of S0 2 and sulfate. Therefore, a priority for future research should be 
to reduce the large uncertainties associated with the wet physical and chemical processes. 
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Fli*’ important roles of sulfate aerosol in climate 
change. atmospheric chemistry, and environmental 
health have been well recognized in recent years. Sul- 
fate aerosol is one of the major aerosol types in the 
troposphere with a dominant anthropogenic compo- 
nent. It affects the Earth's radiative balance directly 
by scattering solar radiation and indirectly by form- 
ing new clouds and modifying cloud properties. It 
also provides surfaces for heterogeneous reactions to 
take place, thus altering the concentrations of many 
important atmospheric species. Sulfate can also inter- 
act with other types of aerosols, such as dust and car- 
bonaceous aerosols, to modify their hygroscopic prop- 
erties when internally mixed with them. The funda- 
mental step toward quantifying all the direct and indi- 
rect effects of sulfate aerosol is determining its spatial 
and temporal distributions and the various processes 
that control the distributions. 

There have been numerous observational data of 
aerosols and their precursors obtained at ground sites, 
in field campaigns, and from satellite measurements. 
However, measurements at the surface or in field cam- 
paigns are limited in spatial or temporal coverage, 
while satellite observations are limited in measurable 
quantities. Therefore, a global model is needed to in- 
tegrate the space-borne, air-borne, and ground-based 
data in order to interpret the data in a broader con- 
text. In fact, several global models have been used to 
study the tropospheric sulfur cycle since 1991 [ e.g ., 
Langner and Rodhe, , 1991; Pham et al , 1995; Fe- 
ichter et al , 1996; Chin et al , 1996; Chuang et al , 
1997; Roelofs et al } 1998; Koch et al, 1999; Barth 
et al, 2000; Rasch et al, 2000]. Almost all the 
published global sulfur models were either driven bv 
the off-line meteorological fields generated in general 
circulation models (GCM), or were coupled on-line 
with the GCM. Although these model studies have 
helped to ad\ance our understanding of tropospheric 
sulfur cycle, it is often difficult for them to explain 
the observed day-to-day and year-to-year variability, 
let alone to interpret in-situ data from fields cam- 
paigns. This is mainly because the results from the 
GCM models in general represent multi-year values 
a\ fudged over a large area, which is inappropriate for 
comparisons with observations in a specific time. 

Here we introduce the Georgia Tech/Goddard Global 
Ozone Chemistry Aerosol Radiation and Transport 
( OOC ART ) model, which can be potentially the most 
suitable um) 1 to link the satellite* and in-situ obser- 


vations. The main advantage of the model, which 
is also the main difference between this model and 
the previously published models, is that the GO- 
CART model is driven by the assimilated meteorolog- 
ical fields, which are generated in the Goddard Earth 
Observing System Data Assimilation System (GEOS 
DAS). This type of model is therefore appropriate 
for interpreting measurements for a specific period of 
time. And because the GOCART model is a global 
scale model, it is also convenient to use in analyzing 
satellite data and conducting global assessments. 

In this paper, we provide a detailed description of 
the model components used for simulating the tropo- 
spheric sulfur cycle (section 2). The global distribu- 
tions and 6-year budgets for sulfate and its precursors 
are presented (section 3), and the anthropogenic con- 
tribution to the sulfate burden is discussed (section 
5). Results from our model are compared with those 
from two most recent model studies [Koch et al, 1999; 
Barth et al, 2000 and Rasch et al, 2000] (section 4 
and 5). A detailed evaluation of the model results 
with observations and budgets for several continen- 
tal and oceanic regions are presented in a companion 
paper [Chin et al, this issue]. It is noted that in addi- 
tion to sulfate, other aerosol components are also sim- 
ulated in the GOCART model, which include dust ( P 
Ginoux et al, manuscript in preparation. 2000). car- 
bonaceous. and sea salt aerosols (in progress). With 
all the major aerosols simulated, we will be able to 
compare the aerosol properties generated in the model 
with those retrieved from the satellite measurements, 
and apply the model to global aerosol analysis. 

2. Model Description 

2.1. Model Framework 

The GOCART model uses the GEOS DAS as- 
similated meteorological data [Schubert et al. 1993]. 
The spatial resolution of the model is the same as 
in the GEOS DAS. which has a horizontal resolution 
of 2 latitude b\ 2. 5° longitude. The vertical resolu- 
tion varies with different versions of the GEOS DAS. 
There are 20 vertical sigma levels in version 1 (GEOS- 
1. available for the period of January 1980-November 
1995t. extending from the surface to 10 mb [ Allen 
et al . 1996: Chin et al. 1998] . In version GEOS- 
1.3 [available from April 1995 to November 1997). 
designed to support the STRAT (Stratospheric Trac- 
ers ot Atmospheric Transport; mission, there are 46 
vertical levels with approximately 26 of them in the 
stratosphere and the model top at 0.1 mb. In our tro- 
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pospheric simulation. we have aggregated the top 23 
levels (from 40 mb to 0.1 mb) in the GEOS-1.3 to 3 
levels and kept the lowest 23 levels (from surface to 40 
mb) as the same resolution as in the GEOS-1.3 such 
that the total number of model vertical levels is 26. 
The lowest 5 layers in both GEOS-1 and GEOS-1.3 
are centered at approximately 50, 250. 600, 1100, and 
1S00 meters above the surface. Newer versions of the 
GEOS DAS data, for example, GEOS-2 and GEOS- 
3, with higher vertical or horizontal resolutions have 
become available for the time periods after November 
1997. 

The GEOS DAS meteorological data contain not 
only prognostic fields, such as horizontal winds, tem- 
perature, and pressure, but also extensive diagnostic 
fields, such as cloud mass flux, surface precipitation 
rates, boundary layer height, and surface roughness. 
Table 1 lists the GEOS DAS archived prognostic and 
diagnostic fields used in our sulfur simulations. 

We present in this paper a 6-year simulation from 
1989 to 1994. Four sulfur species are simulated in 
the model: dimethylsulfide (DMS), S0 2 , sulfate, and 
methanesulfonic acid (MSA). There are seven mod- 
ules representing atmospheric processes of these sul- 
fur species: emission, chemistry, advection, convec- 
tion, diffusion, dry deposition, and wet deposition. 
The model solves the continuity equation using the 
method of operator splitting. The model time step 
is 20 minutes for advection, convection, and diffu- 
sion, and 60 minutes for the other processes. The 
instantaneous meteorological fields in Table 1 are lin- 
early interpolated to the model time. Initialization 
was done for the last 3 months of 1988, starting from 
low concentrations (0.1 ppt) for all four sulfur species. 

2.2. Transport 

The advection and convection schemes used in 
the model have been described in detail elsewhere 
[Allen et ai, 1996]. Here, briefly, advection is com- 
puted by a flux-form semi-Lagrangian method [Lin 
and Rood 1996]. Moist convection is parameterized 
using archived cloud mass flux fields from the GEOS 
DAS. In the previous model studies using the GEOS 
DAS fields, the boundary layer mixing was param- 
eterized such that a fixed fraction of material was 
uniformly mixed within the boundary layer [Allen et 
ui r 1996: Chin et ai . L998j. It was found very diffi- 
cult to choose a universal value of the mixing fraction 
since it does not reflect the boundary layer turbulence 
[Chin et al., 1998] . In the COCART model, the 
boundary layer turbulent mixing is computed using 


a second-order closure scheme [Helfand and Labraga , 
1988], which was also used in the GEOS DAS analysis 
for heat and moisture turbulent mixing [Takacs et ai, 
1994]. The scheme takes into account both growing 
and decaying turbulence. The turbulent diffusion co- 
efficient is a function of the turbulent kinetic energy, 
the buoyancy and wind shear parameters. 


2.3. Sulfur Emissions 


The GOCART model includes emissions of DMS 
from the ocean, S0 2 and sulfate from anthropogenic 
activities, and S0 2 from biomass burning, aircraft ex- 
haust, and volcanic eruptions. Figure 1 shows an an- 
nually averaged emission flux from anthropogenic and 
natural sources (DMS and volcanic S0 2 ) for 1990. 


Figur 


Anthropogenic emissions are taken from the Emis- 
sion Database for Global Atmospheric Research (EDGAR) 
for the year of 1990 [Olivier et al . , 1996]. The annual 
total emission rate is 72.8 Tg S yr” 1 which includes 
emissions from industrial processes (59.3 Tg S yr x ), 
residential and commercial consumptions (8.5 Tg S 
yr” 1 ), and transportation (road, rail, and shipping, 

5.0 Tg S yr -1 ). The fraction of direct sulfate emis- 
sion has been estimated from 1.4% to 5% of the total 
emission [Benkovitz et ai, 1996]; we assume here a 
fraction of 5% for Europe and 3% for elsewhere. The 
rest is emitted as S0 2 . Emission rates are assumed to 
be constant throughout the year except for Europe, 
where a seasonal variation is imposed such that the 
emission rates are maximum in winter (30% higher 
than the annual average) and minimum in summer 
(30% lower than the annual average). This seasonal 
variation reflects mainly the demand for domestic 
heating [Sandnes and Styve , 1992]. 

Emission of DMS from the ocean is calculated as 
a product of the seawater DMS concentration and 
sea- to- air transfer velocity. Monthly averaged surface 
seawater DMS concentrations in l 3 x V grid resolu- 
tion are taken from Kettle et ai [1999]. This sea- 
water DMS concentration map is generated based on 
the compilation of a database of over 15,000 mea- 
surements around the globe. The transfer velocity of 
DMS is computed using an empirical formula from 
Liss and Merlivat [1986]. which assumes linear rela- 
tionships between the transfer velocity and the IO- 
meter wind speed. Diffusion of DMS within the ocean 
surface water is taken into account as a function of sea 
surface temperature [Saltzman et al. , 1993]. The 10-m 
winds used in the model are the remote sensing data 
from the Special Sensor Microwave Imager (SSM/I) 
operated on a series of satellites in the Defense Mete- 
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om logical Satellite Program [Ai/a.s et ai, 1996]. The 
SSM/I winds have been found to represent accurately 
the local observations [Chin et al., 1993]. It has been 
noted that there could be a factor of 2 or more differ- 
ences in the transfer velocity calculated from different 
formulae [e.g., Srnethie et al. f 1985; Wannmkhof et 
al. 1992; Erickson et al. 1993], and a single param- 
eterization of transfer velocity based on wind speed 
alone is not sufficient to describe DMS flux from the 
different regions of the oceans [Chin et ai, 1998]. 

Volcanic sources of S0 2 include emissions from 
both continuously active and sporadically erupting 
volcanoes. The continuous volcanic emissions are 
taken from a database of Andres and Kasgnoc [1998]. 
The database includes S0 2 released from 49 volcanoes 
which have been continuously active over the last 25 
years with an emission rate of 4.8 Tg S yr~ ! . We as- 
sume that S0 2 is injected at a constant rate within 
1 km above the crater altitudes. For the sporadically 
erupting volcanoes, we use the volcanic database from 
the Smithsonian Global Volcanism Program [Simkin 
and Siebert, 1994] which has documented the loca- 
tions. erupting dates and duration, and the volcanic 
explosivity index (YEI) up to 1994. We then use the 
\ El to estimate the volcanic cloud height [Simkin and 
Siebert , 1994], and obtain the amount of S0 2 emit- 
ted to the atmosphere by a relationship between VEI 
and SO-2 flux [ Schnetzler et ai. 1997]. When they 
become available, satellite observed volcanic emission 
data from the Total Ozone Monitoring Spectroscopy 
(TOMS) instrument [Bluth et al . , 1997] are used to 
replace the calculated emission rates. We further 
assume that S0 2 is injected within a slab which is 
located at the top portion of the erupting volcanic 
cloud with a thickness of 1/3 of the cloud column [L. 
Glaze, personal communication, 1998]. This assump- 
tion is based on the observations of plume height and 
thickness after eruption [e.g., McCormick et al . . 1995] 
and the results from volcanic plume dispersion models 
[e.g.. Suzuki, 1983]. 

Other sources of S0 2 in the model include biomass 
burning (2.3 Tg S yr" 1 ) and aircraft emissions (0.07 
Tg S yr ). Seasonal biomass burning emissions are 
from Spiro et al. [1992]. Aircraft emission is calcu- 
lated based on the monthly averaged fuel consump- 
tion inventory for 1992 from NASA's Atmospheric Ef- 
fect.'! of Aviation Project (AEAP), assuming an emis- 
sion index of 1.0. i.e., 1 g SO_> emitted per kg fuel 
burned [IVeisenstein et al.. 19961. 

In our 6-year simulation presented in this paper, 
we have* used the same seasonal emissions from an- 


thropogenic, biomass burning, aircraft, and contin- 
uously active volcanic sources for every year. The 
only interannuallv variable sources are the emissions 
from sporadically erupting volcanoes (based on docu- 
mented events) and DMS from the ocean (due to the 
change of surface wind speeds). 


2.4. Chemistry 


Chemical reactions included in the model are: DMS 
oxidation by OH during the day to form S0 2 and 
MSA, and by nitrate radicals (N0 3 ) at night to’ form 
S0 2 ; S0 2 oxidation by OH in air and by H 2 0 2 in 
cloud to form sulfate. Reaction rates are taken from 
DeMore et al. [1997], The yields of S0 2 and MSA 
from DMS+OH reaction are assumed to be the same 
as in Chin et al. [1996]. i.e.. 100% S0 2 from the ab- 
straction channel, and 75% S0 2 and 25% MSA from 
the addition channel. We prescribe concentrations 
of OH, XO3, and H 2 0 2 from the monthly averaged 
fields generated in the IMAGES model [Muller and 
Brasseur, 1995]. Figure 2 plots the zonallv averaged 
concentrations of OH and H 2 0 2 for January and July. 
A diurnal variation of OH concentrations is imposed 
bj scaling the average OH fields to the cosine of solar 
zenith angle. Since the concentrations of N0 3 over 
the ocean at night are always orders of magnitude 
higher than those during the day. they are assumed 
to be zero in the daytime and are evenly distributed 
over the night. 

Because cloud water content is not available in 
GEOS-1 and GEOS 1.3, we parameterize the in-cloud 
oxidation of S0 2 by H 2 0 2 as a function of cloud frac- 
tion, following Chin et al. [1996], Cloud fraction f e 
for each grid box is assumed as an empirical function 
of the relative humidity in that grid box. following 
Sundqvist et al. [1989]: 
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where r is the relative humidity ami r 0 is the thresh- 
old relative humidity for condensation specified as a 
function of pressure [.Yu and Krueger, 1991], Within 
the cloud fraction we assume that the formation of 
sulfate is determined by the concentration of the lim- 
iting reagent, i.e., the lesser amount between SO-, 
and H 2 0 2 . During the chemistry time step f 1-hour) 
H-’O- is depleted as the aqueous phase reaction of 
S0 2 -rH 2 0 2 can be completed in less than one hour 
[Damn et ai. 198-4]. The recovery time of H 2 0 2 varies 
considerably in the literature, from instantaneous re- 
plenishment to a day in winter [Koch et ai. 1999}. 
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Hero we assume than Hj()_» is regenerated to its pre- 
scribed value every 3 hours, similar to the time scale 
used in Chin et al. [1996]. 

2.5. Dry Deposition 

Dry deposition velocities for S0-2 ? sulfate, and 
MSA are calculated using the resistance-in-series scheme 
[Wesely and Hicks , 1977]. In this scheme, dry depo- 
sition velocities are determined as a reciprocal of the 
sum of aerodynamic resistance, sub-layer resistance, 
and surface resistance. The aerodynamic resistance is 
taken from the GEOS DAS archive, which is a prod- 
uct of the exchange coefficient for heat and moisture 
and the surface friction velocity. The sub-layer and 
surface resistance for SO? and sulfate are calculated 
using the formulation of Walcek et al. [1986] and We- 
sely [1989]. The dry deposition velocity of MSA is 
assumed to be the same as that of sulfate. We im- 
pose a minimum SO 2 dry deposition velocity of 0.2 
cm s' 1 over the ice and snow and in the polar regions 
[Voldner et al. , 1986; Tarrason and Iversen , 1998]. 
Typically, the diurnally averaged dry deposition ve- 
locity for S0 2 over the land is 0.2-0. 4 cm s"\ but it 
varies significantly over the ocean, from 0.6 to 0.8 cm 
s" 1 under stable conditions to 1 to 2 cm s' 1 under 
unstable conditions. For sulfate, a value of 0.08-0.12 
cm s“ l is found over the oceans and 0. 1-0.3 over the 
land except at latitudes higher than 60° in the win- 
ter season (0.01-0.05 cm s' 1 ). These values are in 
general consistent with the data from limited direct 
measurements and with other calculated values [e.g., 
Voldner et al, 1986; Walcek et al, 1986; Ganzeveld 
et al. 1998 and references therein]. 

2.6. Wet Scavenging 

Wet scavenging of soluble species in the model in- 
cludes rainout (in-cloud precipitation) and washout 
(below cloud precipitation) in large-scale precipita- 
tion and in deep convective cloud updraft. The 
GEOS DAS diagnoses the total precipitation at the 
ground as an column integral of specific humidity 
change due to moist processes [Takacs et al. 1994]. 
Here we normalize the precipitation rates from the 
GEOS DAS to those from an observation based data 
product, which is a merged dataset combining satel- 
lite observations, ground station rain gauge measure- 
ments, and the GEOS DAS precipitation fields [P. 
Houser, manuscript in preparation. 2000: Huffman et 
al. 1997]. Distribution of large-scale precipitation in 
a vertical column is estimated based on the specific 


humidity changes diagnosed in the GEOS DAS: 

PiJk) = cMt(k) 

where Pu(k) is the large-scale precipitation rate at 
level k, Q u and Q t are respectively the column inte- 
grated specific humidity change due to large-scale or 
total (large-scale and convective) moist process, c is 
the ratio of the precipitation rate in the merged prod- 
uct to that in the GEOS DAS, and i\qt{k) is the total 
specific humidity change at level k where a negative 
value indicates a precipitation and a positive value 
implies an evaporation. 

Removal of sulfate and MSA by large-scale rain is 
calculated as a first-order loss process using parame- 
ters of Giorgi and Chameides [1986]. The change of 
aerosol mixing ratio within a model time step is: 

Ax(fc) = X (A - 1) 

where \{k) is the mixing ratio at level k, f(k) is the 
fraction of the grid box experiencing precipitation, 
8(k) is the frequency of cloud to rain conversion, and 
At is the duration of precipitation, which is equal to 
the wet scavenging time step for large-scale rain. The 
values of f(k) and 3{k) are defined by the precipita- 
tion amount at each grid box and by a typical liquid 
water content for large-scale precipitation [Giorgi and 
Chameides , 1986]. 

Washout between the cloud layers or below the 
lowest cloud level is also computed as a first-order 
loss process, similar to the treatment of rainout. In 
this case, the fraction of a grid box with precipita- 
tion is determined by the largest value of / from the 
overhead rainy grid box. and 3 is assumed to be 0.1 
mm' 1 normalized to the precipitation rate [Dana and 
Hales , 1976]. A fraction of soluble species between or 
below clouds releases into the grid box if evaporation 
(A^ > 0) occurs. This fraction is assumed to be the 
same as that of evaporated water. 

It has been found in previous model investigations 
as well as in field studies that soluble species are scav- 
enged efficiently within the convective cloud updraft 
[Balkanski et al. 1993; Cohan et al. , 1999]. Adapting 
the principle of Balkanski et al [1993]. wo couple the 
convective scavenging with the moist convection pro- 
cess in our model, and use a scavenging efficiency of 
0.4 km -1 for soluble aerosol species. 

We use the same method for SO 2 wet scavenging 
as that described in Chin et al. [1996]: we define a 
soluble fraction of SOj as limited by the availability of 
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HjO > in tiie precipitating grid box. and scavenge the 
soluble SO, at the same rate as sulfate. When evap- 
oration occurs, a fraction of dissolved SO, returns to 
the grid box as sulfate. 

3. Global Budget and Distributions 

3.1. Summary of Global Budget 

Summary of a 6-year budget of 1989-1994 is pre- 
sented in Figure 3. Before we discuss the budget, 
rve shall clarify the terms used in our wet removal 
and aqueous-phase oxidation budgets, since they can 
sometimes cause confusion. Here, the term “wet scav- 
enging refers to the loss of a particular tracer in the 
wet process described in section 2.6 regardless of its 
transformation within the rain water. With that in 
mind, we count the amount of SO, scavenged and 
subsequently converted to sulfate in the rain water as 
a term of wet scavenging of SO,, not sulfate. And we 
do not record this amount as a part of “in-cloud sul- 
fate production (except for the fraction returned to 
the atmosphere during the evaporation of raindrops), 
because the production of sulfate from the dissolved 
SOj in rainwater does not contribute to either the 
sulfate burden in the atmosphere or the removal of 
sulfate from the atmosphere. While it seems just a 
labeling issue for the sulfur budget, counting the wet 
scavenging of SO, as a loss of sulfate can lead to an 
underestimation of atmospheric sulfate lifetime, since 
the lifetime is simply the ratio of atmospheric burden 
to the loss rate. 

_As shown m Figure 3, the anthropogenic emission 
is ( 5 Tg S yr l , which includes emissions from indus- 
trial activities, fuel combustion, ship, and aircraft, 
as well as from biomass burning. Biogenic emission 
of DMS from the ocean varies from 13.3 to 15.0 Tg 
S }i , reflecting the changes in the surface wind 
speeds. Volcanic emissions are also fairlv constant 
from year to year (5.4-6.0 Tg S yr" 1 ), except 1991 
when a major volcanic eruption of Mt. Pinatubo oc- 
curred in June, injecting about 10 Tg S (or 20 Mton 
SO,) into the atmosphere. Total volcanic emission for 
1991 is 19.6 Tg S yr -1 . In a normal year (e.g., with- 
out major volcanic eruptions), the fraction of sulfur 
emitted from natural sources (biogenic and volcanic) 
is about 20% of the total emission of 94 Tg S yr -1 . 

In-cloud oxidation of SO, is responsible for about 
64 T. of total sulfate production in a normal year, while 
m-air oxidation accounts for the rest 36 %. In con- 
ti.ist. less than half of the sulfate production in 1991 
takes place in-cloud, because the Pinatubo eruption 


injects most SO, into the stratosphere where the gas- 
phase reaction with OH is the onlv mechanism in 
the model to convert SO, to sulfate. Dry deposition 
and wet scavenging remove roughly the same amount 
ot sulfur from the atmosphere (45-55 Tg S yr -1 ). 
W hile dry deposition is the most important loss of 
SO, (45%) followed by in-cloud oxidation (27%), wet 
scavenging eliminates 90% of sulfate produced in the 
atmosphere. The lifetime is 1.8 days for SO, and 5.8 
days for sulfate in a normal vear. 


u.uajcu aunuspneric ourcien tor SO, 

is 0.42-0.48 Tg S except 1991. The SO, burden is 1.6 
Tg S in 1991, with most of it residing in the strato- 
sphere, i.e., above 100-120 mb in the model (Fig- 
ure 3). While S0 2 returns to its normal level rather 
quickly after the Pinatubo eruption (e-folding time 
about 1 month), it takes much longer for sulfate to 
relax back to its normal level. As illustrated in Fig- 
ure 3. three years after the Pinatubo eruption, total 
sulfate burden in 1994 (0.98 Tg S) is still significantly 
higher than its pre-Pinatubo value in 1989-1990 (0 63 
T S S). 


The only removal process for DMS is its oxida- 
tion in the atmosphere. Globally, nearly 90% of DMS 
emitted from the ocean is oxidized by OH during the 
day: only 10% is lost at night via reaction with N0 3 . 
The stable products from DMS oxidation are 89% 
S0 2 . which can be further oxidized to sulfate, and 
11% MSA, which is removed by wet (91%) and dry 
(9%) depositions. The atmospheric burden for DMS 
is 0.0/2-0.080 Tg S, and that for MSA is 0.027-0.032 
Tg S. The lifetime is 1.9-2.2 days for DMS and 6.S-7.2 
days for MSA. 


We define a term of the sulfate production effi- 
ciency as the amount of sulfate produced relative to 
the total amount of SO, emitted and produced in 
the atmosphere. The sulfate production efficiency is 
a direct measure of the effectiveness of SO, oxida- 
tion versus the dry and wet removal of SO,. We find 
in our model a typical production efficiency value of 
0.41-0.42, which indicates that only less than half of 
the SO, contributes to sulfate production in the atmo- 
sphere, and the rest is mainly deposited to the surface. 
In the Pinatubo eruption year of 1991. however, the 
sulfate production efficiency increases to 0.49. reflect- 
ing that the SO, released at high altitudes produces 

sulfate much more effectively than that emitted near 
tiie surface. 



3.2. Global Distributions 

To present some general features simulated in the 
model, we plot in Figure 4 global distributions of S0 2 , 
sulfate, DMS. and MSA at the surface (Figure 4a) and 
at 500 mb (Figure 4b) for the pre-Pinatubo year of 
1990. Concentrations shown in Figure 4 are average 
values for 2 seasons: December, January. February 
(DJF), and June, July, August (JJA). High surface 
concentrations of S0 2 and sulfate are found in regions 
of high anthropogenic emissions for both seasons, as 
expected. The major contrast between DJF and JJA 
is the strong advection of pollutants from the mid- 
latitude source regions to the Arctic circle in DJF. 
While S0 2 concentrations are higher in the winter 
than in the summer, the reverse is true for sulfate, 
because of the seasonal variation of S0 2 oxidation 
rates. Globally, the sulfate production efficiency in 
January is only 0.27 whereas in July it is 0.48 in 1990. 

The distribution of DMS at the surface closely re- 
sembles that in seawater. Very high surface air con- 
centrations of DMS (500-2000 ppt) are produced in 
the model near 60° latitude in the summer hemi- 
sphere. These elevated concentrations are directly 
related to the high DMS emission flux (10-50 ^mol 
m" J day" 1 ), a product of very high seawater DMS 
concentrations and strong surface winds. \V hile the 
surface DMS concentrations at these latitudes in sum- 
mer seem too high compared with some measure- 
ments near the Antarctic coast (< 800 ppt. Staubes 
and Georgii , 1993: Berresheim et al , 1998), the 
model calculated concentrations of the DMS oxida- 
tion products, sulfate and MSA, at high latitude sites 
(e.g.. Palmer Station and Mawson in Antarctica, and 
Haemey in Iceland) agree with the observations to 
within 40% [Chin et al, this issue]. This apparent in- 
consistency needs to be further investigated. Finally, 
as expected. MSA surface distribution is similar to 
that of DMS. 

One common feature for all sulfur species at 500 
mb (Figure 4b) is that they are better mixed zon- 
allv and the concentrations are 1 to 2 magnitudes 
lower than that at the surface, reflecting their rela- 
tively short lifetimes (several days). DMS concentra- 
tions at 500 mb are higher in the winter hemisphere 
than in the summer hemisphere, opposite to the pat- 
tern found at the surface. This is because the slow 
oxidation rates of DMS in winter allow DMS to be 
transported to higher altitude's. On the other hand, 
the DMS loss rates are higher in summer than in win- 
ter due to the higher OH concentrations: thus only 
a small fraction of DMS escape's from the boundary 


layer in summer despite the higher emission rates. 

The annually averaged zonal mean distributions 
of S0 2 . sulfate, DMS, and MSA for 1990 are shown 
in Figure 5. As expected, both SO_> and sulfate ex- 
hibit high concentrations in the northern hemisphere. 
In the tropics. DMS is pumped to the upper tropo- 
sphere by the deep convective process. Interestingly, 
the same process is also responsible for the low sul- 
fate concentration in the middle to upper troposphere 
over the tropics, thanks to the efficient wet scaveng- 
ing of sulfate in cloud convection. This feature also 
appeared in the model simulations of Feichtev et al. 
[1996] and Koch et al [1999], but was lacking in some 
other models (e.g., Chin et al, 1996; Barth et al, 
1999), depending on the convective process in the me- 
teorological data and the efficiency of in-cloud scav- 
enging in different models. 

We plot in Figure 6 the column total sulfate sources 
and sinks in the 1990 simulation as a function of lat- 
itude. It can be seen that in-cloud oxidation of S0 2 
is the most important source of sulfate, especially 
at high latitudes (60 C N and higher in the northern 
hemisphere, and between 40°S and 70°S in the south- 
ern hemisphere), where in-cloud oxidation contributes 
80-90% of the total sulfate source. While scaveng- 
ing by large-scale rain dominates the sulfate loss over 
mid- and high latitudes, wet convective scavenging 
dominates over the tropics and subtropics. Dry depo- 
sition in general accounts for less than 20% of sulfate 
loss at all latitudes except in the polar regions. 


Figu: 


Figui 


4. Comparisons with Other Global 
Model Studies 


We focus here on the comparison of sulfur budget 
in the GOCART model simulation for 1990 with two 
of the most recent global model studies: Koch et al. 
11999]. using the Goddard Institute for Space Studies' 
general circulation model version II (GISS GCM), 
and Barth et al [2000] and Rasch et al [2000 j, us- 
ing the NCAR Community Climate Model (NCAR 
CCM3). Intercomparisons involving other earlier 
models (e.g., Langner and Rodhe. 1991: Pham et al, 
1995: Feichtev et al , 1996: Chin et al , 1996) have 
been presented in previous model studies [ Chin et al, 
1996: Koch et al. 1999: Rasch et al. 2000] and will 
not be discussed in derail here. We will only summa- 
rize tin' major differences between this study and an 
earlier work [Chin et al.. 1996]. 

Table 2 summarizes the comparison of sulfur bud- 
gets among the GOCART. GISS. and NCAR mod- 
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'■I s - <),|r total emission (93.9 Tg S vi- 1 ) is higher 

rhan thiU both GISS and NCAR models (83 Tg S 
1 ~ 1 - • Tll 's is because we use the EDGAR database 
ot 1990 emission inventory (72.3 Tg S yr" 1 ) which also 
includes emissions from shipping and landuse. while 
the GISS and NCAR models use the GEIA emission 
inventory (67 Tg S yr" 1 ) for the 1985 emission sce- 
nario [Benkovitz et d.. 1996], The biomass burning 
emission in our model is the same as that in the GISS 
model (2d3 Tg S yr -1 ), but the volcanic emission (5.5 
Tg S yr“i) is much higher than that in the GISS 
model (3.5 Tg S yr -1 ), because we include emissions 
from both continuously active and sporadically erupt- 
ing volcanoes, while Koch et d. [1999] considers only 
non-eruptive volcanoes. Volcanic and biomass burn- 
ing emissions are not included in the NCAR model. 
Emission of DMS calculated in our model is 13.3 Tg 
S yr- 1 for 1990, about 30% higher than that in the 
GISS model (10.7 Tg S yr -1 ), even though both mo- 
dles use the same formula for transfer velocity and 
the same DMS seawater concentrations in calculat- 
ing DMS emission rates. This difference may be at- 
tributed to the lower 10-meter wind speeds in the 
GISS GCM. The NCAR model calculates a total DMS 
emission of 15.5 Tg S yr -1 , based on the latitudinal 
bands of DMS flux in Bates et al. [1992] and the 
distribution of the ocean color in the remote sensing 
products. 

Our model estimates an equal amount of sulfur be- 
ing removed by dry deposition and wet scavenging 
(50% for each process), while a slightlv higher frac- 
tion of dry deposition (54% dry, 46% wet) is obtained 
in the GISS model and the wet removal is about twice 
as effective as dry deposition in the NCAR model. 

Total sulfate production from S0 2 oxidation is 38.5 
Tg S yr 1 in our model. Although this value is the 
lowest among the three models, remember that we do 
not count the S0 2 loss in wet scavenging as a part of 
sulfate production while both the GISS and NCAR 
models do. As we stated in the previous section, even 
though the amount of S0 2 scavenged by the rain is 
subsequently converted to sulfate in rainwater, this 
process does not play a role in determining the at- 
mospheric sulfate concentration or removal: we thus 
consider the in-rain sulfate production as ineffective. 
Should we include the S0 2 wet scavenging as a part 
of sulfate production, the value would be 49.1 Tg S 
yr' 1 which is between the GISS and NCAR model. 

As w<* have shown in the previous section. S9% 
of the DMS emitted from the ocean produces SO_> 
(11.9 Tg S yr *). Of this amount, 87% is produced 


via DMS-fOII and 13% via DMS + NO.j. The SO-> 
production is more efficient in the GISS model (93%)’ 
whereas the NCAR model assumes S0 2 as the onlv 
DMS oxidation product. 

The GISS model has the highest SO_> and sulfate 
burden among the three models, which was attributed 
to the use of prognostic H 2 0 2 and an insufficient en- 
trainment of H 2 0 2 from the cloud base to oxidize SO ? 
in highly polluted regions [ Koch et d.. 1999], On the 
other hand, the NCAR model also uses prognostic 
H->°2 but shows the lowest S0 2 and sulfate burden. 
The cause of the discrepancy is likely a combination 
of the differences in cloud processing, oxidant concen- 
trations, precipitation rates, among others, between 
the models. The DMS burden and lifetime in our 
model is 22% and 41% higher than those in the NCAR 
model, even though both models have used the same 
prescribed OH and N0 3 fields for DMS oxidation. 
Possible explanations include the difference in DMS 
emission rates, which are higher in our model at high 
latitudes where DMS oxidation is slower The life- 
time of DMS in the GISS model is verv close to that 
in our model, although the DMS burden is lower in 

the GISS model, probably due to the lower emission 
rates. 

The lifetime of 2.6 days for S0 2 in the GISS model 
is about 40% longer than the ones in both our model 
(1.8 days) and the NCAR model (1.9 davs), reflect- 
ing a slower removal rate in the GISS model. Re- 
garding the lifetime of sulfate, we estimate a value 
of 5.8 days with respect to the sulfate drv deposition 
and wet scavenging. If we had included the amount 
of SO> lost by wet scavenging as a sink of sulfate, 
as the GISS and NCAR models have, then the sul- 
fate lifetime in our model would be 4.6 days, which 
is lower than the GISS model but higher than the 
NCAR model. The point we try to make here is that 
the atmospheric sulfate residence time is underesti- 
mated when wet scavenging of S0 2 is included as a 
loss term of sulfate. The lifetime of DMS is similar to 
that in the GISS model but higher than that in the 
NCAR model, whereas the lifetime of MSA is slRhtlv 
lower than that in the GISS model (MSA is not°sim- 
uiatecl in the XCAR model). 

To examine the differences between the models in 
loss rates for individual sinks, we list in Table 2 the 
loss frequencies for each process, defined as the SO, 
or sulfate atmospheric burden divided by their indi- 
vidual removal rates. While dry removal processes 
for SO, (dry deposition and in-air oxidation) are the 
most efficient in our model, wet processes (aqueous 
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phase product, ion and wet removal) are the most ef- 
fective in the N CAR model but the least effective in 
the GISS model. The effectiveness of the wet process 
is almost inversely proportional to the SO_> and sul- 
fate burden, which are the lowest in the NCAR model 
and highest in the GISS model. For example, a ratio 
of the sulfate burden between the GOCART, GISS, 
and NCAR models are 1:1. 3:0. 9. that of the SOo bur- 
den is 1:1. 2:0. 9, which is also the ratio of the inverse 
of the total sulfur wet deposition rate. Interestingly, 
the sulfate burden is also inversely proportional to the 
wet production rate of sulfate (in-cloud and in-rain); 
a ratio of 1:1. 1:0.9 is found between the three models. 
These linear relations clearly confirm the importance 
of the wet processes in determining the sulfur bur- 
den in the atmosphere. Therefore, emphasis should 
be given to improving the wet physical and chemi- 
cal processes and validating the parameters used in 
modeling these processes, such as cloud distribution, 
cloud fractions, precipitation amount, scavenging ef- 
ficiencies, and the aqueous phase oxidation rates. 

Comparing our zonally averaged concentrations of 
SO>. sulfate, and DMS in Figure 5 with those re- 
ported in Koch et ai [1999] and Barth et ai [2000], 
we find that the SO 2 concentrations in the GISS 
model are significantly higher than that in both the 
NCAR and our model in the lower troposphere. For 
example, a 500 ppt SO 2 contour in the GISS model 
reached 700 mb and extended from 27°N to 75°N, 
while this contour line is confined below* 850 mb and 
at the latitudes between 25 s N and 60° N in both the 
NCAR and our models. A similar difference in sulfate 
distribution is also found. Our extratropical DMS 
zonal distribution resembles that in the GISS model 
w*ith a symmetric distribution between the northern 
and southern hemispheres. However, in the tropi- 
cal upper troposphere, there is a second maximum of 
DMS in our model with a concentration of 5-10 ppt, 
a feature which is very similar to the NCAR model 
but not obvious in the GISS model. 

The results from the GOCART model differ from 
that in Chin et al [1996] (using the Harvard/GISS 
GCM II model) in a number of ways. The ma- 
jor difference is in the sulfate vertical distributions. 
The zonally averaged sulfate distribution in Figure 5 
shows a much less vertical gradient than that in Chin 
et ai [1996]. The difference is attributed mainly 
to a much more efficient wet scavenging in Chin et 
al. [ 1996] . partly due to the excessive wet convection 
over some regions, and partly due to the higher wet 
scavenging efficiency (100% in deep wet convection). 


The other difference is in DMS oxidation. Chin et al. 
[1996] found that an oxidant in addition to OH and 
NO ; j was needed for DMS oxidation in order to repro- 
duce both DMS and sulfate concentrations observed 
over the remote ocean surface. We do not invoke such 
an oxidant in this study, and our simulated concen- 
trations for all sulfur species are overall consistent 
with the observations over the oceans [Chin et al. , 
this issue]. We attribute this difference to the bet- 
ter parameters in calculating the DMS emission rates 
and the higher (a factor of 2 to 3) OH concentmtions 
over the ocean surface (a factor of 2 to 3) used in this 
study than those used in Chin et ai [1996]. 


5. Anthropogenic Contributions 


We have conducted a model simulation for 1990 
without anthropogenic emissions in order to estimate 
the relative importance of natural and anthropogenic 
sources to the atmospheric sulfate loading. The total 
emission for this case is 18.8 Tg S vr“ l , wffiich in- 
cludes only the DMS and volcanic SO 2 sources. The 
annually averaged column sulfate burden and the an- 
thropogenic contributions for 1990 are shown in Fig- 
ure 7. The anthropogenic fraction of sulfate is more 
than 60% in the northern hemisphere, with more than 
809c over the United States and the Eurasian conti- 
nent. A more widely spread anthropogenic influence 
over the northern hemisphere is found in the GISS 
[Koch et al. 1999] and NCAR [Rasch et ai, 2000] 
models, with more than 80% anthropogenic sulfate 
over the entire area at latitudes north of 10°N. In 
the southern hemisphere, the anthropogenic fraction 
is generally 20-40% over the ocean in our model (Fig- 
ure S). similar to the GISS and NCAR models. 

We find that anthropogenic sources contribute to 
679c of the total sulfate burden in 1990. a fraction 
which is somewhat lower than the anthropogenic sul- 
fur emission fraction of 80%. Figure S shows the 
percentage of zonally averaged anthropogenic contri- 
bution for tw'o seasons, D.JF and JJA, in 1990. As 
can be seen in Figure 3. the anthropogenic sulfate 
dominates the sulfate burdens in the northern hemi- 
sphere but with distinct patterns between D.JF and 
JJA. It disperses horizontally in DJF with the 80% 
contourlino stretched out to the northern polar re- 
gion but confined below 600 mb. By contrast, the 
anthropogenic sulfate is well mixed vertically by the 
frequent convective activities in JJA with the 30% 
contour line extended to the tropopause. Of inter- 
est is that the anthropogenic contribution increases 


Figur 


Figur 
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wirh r.ht? altitude over the mid- to high latitude in the 
southern hemisphere, resulting from the interhemi- 
spheric transport from the northern hemisphere and 
the convective transport from the mid-latitudes. 
When comparing the natural sulfur budget with 
that in the GISS model [ Koch et al, 1999] , a major 
disagreement lies in the lifetimes of S0 2 and sulfate. 
W hile in the GISS model the S0 2 lifetime from a 
natural-source-only run (1.8 days) was shorter than 
that from a full run (2.6 days), we find that the reverse 
is the case in our model: 2.4 days in the natural- 
source-only run and 1.8 days in the full run. As for 
sulfate, the lifetime stayed the same in both natural 
and full simulations in the GISS model, but in our 
model it is longer from the natural simulation (7.2 
days) than that from the full simulation (5.8 days). 
It is expected that S0 2 and sulfate of natural origin 
should have a longer lifetime than the anthropogenic 
ones, because they are not as concentrated near the 
surface, thus not subject to the fast removal by dry 
and wet depositions. 

The anthropogenic contribution to the atmospheric 
sulfate burden from this study, as well as from the 
GISS and NCAR models, is significantly higher than 
that reported by Chin and Jacob [1996]. The lat- 
ter study found that the anthropogenic sources con- 
tributed to only 37% of the sulfate burden, although 
they accounted for 70% of the total sulfur emission. 
This is due to the high sulfate production rates from 
DMS oxidation and more excessive wet scavenging 
near the mid-latitude continents in Chin et al [1996] 
than those in this study. 

6. Conclusions 

We have used the GOCART model to simulate the 
tropospheric sulfur cycle. The model uses the assimi- 
lated meteorological fields from the GEOS DAS. mak- 
ing it potentially the best tool to link the satellite and 
in-situ observations for global analysis. We have in- 
corporated in the model the most updated emission 
inventories of anthropogenic, biogenic, and volcanic 
sources for DMS and SO>. In a typical year with- 
out major volcanic eruptions, we estimate that about 
2u% of the sulfate precursor emission is from natural 
sources (biogenic and volcanic) while 30% is anthro- 
pogenic. In-air and in-cloud oxidation of S0 2 account 
for 36 /c and 64% respectively of the atmospheric sul- 
fate' production. We have estimated a sulfate pro- 
duction efficiency as a ratio of the amount of sulfate 
produced to the total amount of SO_> emitted and pro- 


duced in the atmosphere. A typical production effi- 
ciency value of 0.41-0.42 is found, indicating that gen- 
erally more than half of the S0 2 entering the atmo- 
sphere does not contribute to the sulfate production 
but is either removed by dry deposition or scavenged 
by the rain. We have reported that in 1990 the atmo- 
spheric burdens for S0 2 . sulfate. DMS, and MSA are 
0.43.0.63, 0.073, and 0.023 Tg S, respectivelv, with 
the corresponding lifetimes of 1.8. 5.8. 2.0. and 7.1 
days. 

The anthropogenic contribution to the atmospheric 
sulfate burden is estimated at 67% for 1990, a fraction 
which is somewhat smaller than that of anthropogenic 
emission (80%). While it is horizontally spreading 
out to the northern polar region in DJF, the anthro- 
pogenic contribution is vertically well mixed in JJA, 
with the 80% contour line extended to the tropopause 
over the mid-latitudes in the northern hemisphere. 
We have also shown that major volcanic eruptions can 
significantly change the sulfate formation pathways, 
distributions, abundance, and lifetime. These effects 
are seen in our model simulations from 1989 to 1994, 
a period which includes the major volcanic eruption 
of Mt. Pinatubo in 1991. It has been demonstrated 
that while S0 2 returns to its normal level in only 
a few months after the Pinatubo eruption, it takes 
several years for sulfate to relax back to its normal 
atmospheric loading. 

Our model results of 1990 have been compared 
with two most recent model studies, namely the GISS 
model [ Koch et al . , 1999 ] and the NCAR model 
[Barth et al . , 2000: Rasch et al , 2000]. While the an- 
nual DMS burden in our model is 20-30% larger than 
the other two models, our S0 2 and sulfate burdens are 
lower than those in the GISS model but higher than 
those in the NCAR model. The relative abundance of 
the S0 2 and sulfate burden is almost inversely propor- 
tional to the rate of wet removal and the rate of wet 
production of sulfate. This proportionality shows the 
magnitude of the wet processes in controlling the at- 
mospheric sulfur burden. Therefore, the first priority 
in future research should be to reduce the large uncer- 
tainties associated with the wet physical and chemical 
processes. 
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Figure 1. Annual emissions (mg S nr 2 yr~‘) of 1990 for sulfate precursors from anthropogenic and natural 
(oceanic and volcanic) sources used in the model. 

Figure 2. Zonal averaged distributions of OH and H 2 0 2 for January and Julv from the IMAGES model fA/«// e r 
and Brasseur, 1995]. 1 

Figure 3. Summary of a 6-year sulfur budget (1989-1994) in the GOCART model. Troposphere and stratosphere 
interface at 120-100 mb. 

Figure 4. Distributions of sulfur species (ppt) in for DJF and JJA in 1990 at (a) the surface and (b) 500 mb. 
Figure 5. Zonal distributions of sulfur species (ppt) in the simulation for 1990. 

Figure 6. Sources and sinks for sulfate in the 1990 simulation as a function of latitude. 

Figure 7. Total sulfate column burden (mg S m~ 2 ) and the anthropogenic fraction in the 1990 simulation. 
Figure 8. Zonally averaged anthropogenic sulfate fraction in DJF and JJA for 1990. 



Table 1. GEOS-DAS meteorological fields used in 
GOCART 


GEOS-DAS fields 


Prognostic fields (instantaneous value, 
every 6 hours): 

Surface pressure 
Temperature 
Wind velocity 
Specific humidity 
Surface albedo 

Surface type (land, water, or ice) 

Diagnostic fields (average value, every 
3 or 6 hours): 

Cloud mass flux 
Convective cloud detrainment 
Cloud fraction (column) 

Specific humidity change 

Total precipitation at the surface 

Convective precipitation at the surface 

Aerodynamic resistance 

Surface friction velocity 

Surface roughness length 

Surface air temperature 

Surface sensible heat flux 

Boundary layer depth 

Wind velocity at 10 m 

Turbulent diffusion coefficient 

Xet shortwave radiation at the surface 


a Turbulent diffusion coefficients were not archived in 
the earlier version of GEOS-DAS (before 1997). They 
have been calculated using the archived GEOS DAS fields 
for simulations before 1997. 
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Table 2. Comparison of Sulfur Budget from the GOCART Model with The GISS and NCA.R 
Models. 


Budget Component 

GOCART a 

GISS b 

NCAR C 

Total emission (Tg S yr” 1 ) 

93.9 


83.0 


82.5 


SOo anthropogenic 

70.6 

(75.2%) 

64.6 

(77.8%) 

65.7 

(79.6%) 

SO> biomass burning 

2.3 

(2.4%) 

2.3 

(2.8%) 



SOj volcanic 

5.5 

(5.9%) 

3.5 

(4.2%) 



Sulfate anthropogenic 

2.2 

(2.3%) 

1.9 

(2.3%) 



DMS oceanic 

13.3 

(14.2%) 

10.7 

(12.9%) 

15.5 

(18.8%) 

Total deposition (Tg S yr” 1 ) 

93.0 


83.4 


81.0 


SO2 dry deposition 

41.2 

(44.3%) 

35.5 

(42.6%) 

24.5 

(30.2%) 

SO-> wet scavenging 

10.6 

(11.4%) 

0.2 

(0.2%) 

1.6 

(2.0%) 

Sulfate dry deposition 

5.1 

(5.5%) 

9.2 

(11.0%) 

3.7 

(4.6%) 

Sulfate wet scavenging d 

34.7 

(37.3%) 

37.4 

(44.8%) 

51.2 

(63.2%) 

MSA dry deposition 

0.1 

(0.1%) 

0.2 

(0.2%) 


MSA wet scavenging 

1.3 

(1.4%) 

0.9 

(1.1%) 



SO> production (Tg S yr” 1 ) 

11.9 


10.0 


15.5 


From DMS-f OH 

10.4 

(87.4%) 




From DMS+NO3 

1.5 

(12.6%) 





Sulfate production (Tg S yr" 1 ) 

38.5 


44.7 


53.6 


In-air 

14.0 

(36.4%) 

3.1 

(29.3%) 

9.2 

(17.2%) 

In-cloud d 

24.5 

(63.6%) 

31.6 

(70.7%) 

4.4 

(82.8%) 

Burden (Tg S) 







SO, 

0.43 


0.56 


0.4 


Sulfate 

0.63 


0.73 


0.57 


DMS 

0.073 


0.056 


0.06 


MSA 

0.028 


0.023 



Lifetime (days) 







S0 2 

1.8 


2.6 


1.9 


Sulfate 6 

5.8 

(4.6) 

5.7 


4.0 


DMS 

2.0 

1.9 


1.4 


MSA 

7.1 


7.6 




Loss frequency f (day"’ i ) 







SO2 dry deposition 

0.26 


0.17 


0.17 


SOj in-air oxidation 

0.09 


0.06 


0.06 


SO2 in-cloud oxid.+wet scav. 

0.22 


0.15 


0.31 


Sulfate dry deposition 

0.02 


0.03 


0.02 


Sulfate wet scavenging 

0.15 

(0.2) 

0.14 


0.25 



‘This work, 1900 simulation. 
h Koch et ai, [1999]. 


' Barth et ai , [2000]; Rasch et aL f [2000]. 

d S0 2 wet scavenging is counted as a part of sulfate in-cloud production and sulfate wet scavenging in the 
GISS and NCAR model. See text for details. 

e The numbers in parentheses for the GOCART model are the values that would be if SO 2 wet scavenging 
were considered as a part of sulfate wet deposition term, as treated by the GISS and NCAR models. See 
text for explanation. 

f Loss frequency is defined as the loss rate divided by the burden. 
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